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Abstract 

We study the locahzation of wave functions for one-dimensional Schrodinger Hamihoni- 
ans with random potentials V{x) with short range correlations and large local fluctuations 
such that J dx {V{x)V{0)) = oo. A random supersymmetric Hamiltonian is also considered. 
Depending on how large the fluctuations of V{x) are, we find either new energy dependences 
of the localization length, £ioc oc E/\nE, iioc oc E^^'^ with < /i < 2 or £ioc oc In''"^ E for 
/i > 1, or superlocalization (decay of the wave functions faster than a simple exponential). 

PACS numbers : 72.15.Rn ; 73.20.Fz ; 02.50.-r. 

Introduction.— The phenomenon of Anderson localization [1] in one dimension has been widely 
studied since the pioneering work of Mott & Twose arguing that all states are localized in one- 
dimension (Id) This statement was rigorously proven in Refs. OS]. A general method 
to study the spectral and localization properties of Id random Hamiltonian was proposed in 
Refs. [H [6] : let us consider the one-dimensional Schrodinger Hamiltonian H = —-^ + V{x) 
where V{x) is a random potential with short range correlations. We first study the solution of the 
stationary Schrodinger equation HiJj{x;E) = Eip{x;E) satisfying ip{0;E) = and ijj'{0;E) = 1 
(differenciation with respect to x is denoted by '). We define the Lyapunov exponent (inverse 
localization length ^loc = I/7) as the increase rate [5l[7] 



7(^)= lim ^ / In Ji;{x;E)^ + I; i;' {x; E)A . (1) 
x-tco dx \ \ E I 

Averaging (• • •) is taken over realizations of the random potential. This definition becomes more 
clear if the wave function is parametrized in terms of an envelope and an oscillating part. We 
substitute to the couple of functions (V', V'') the variables (0, ^) according to : 

^l){x\E)= e^^'') sin 61(2;) (2) 
^''(x;^) = A;e«(^) cose(x) , (3) 

where E = k^. We can rewrite the definition of the Lyapunov exponent as jiE) = -^{(.{x))- 
Therefore the Lyapunov exponent gives the rate of the exponential increase of the envelope of 
the wave function. 

At high energy (compared to disorder), oscillations of the wave function occur on the typical 
scale and the Lyapunov exponent is given by [5l [6^] 

-f{E -^oo)c^-^ dx {V{x)V{0)) cos 2kx . (4) 



^ The solution ip{x; E) of the Cauchy problem exists for any value of the energy ; it is used to construct 
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The question of the present article is to discuss the situation where the potential presents 
large fluctuations such that J dx {V{x)V{0)) = oo, which makes (jH inapplicable. The integral 
of the correlation function may diverge for different reasons. First, it may diverge due to long 
range correlations. This corresponds to nonstationary potentials. The localization for self affine 
potentials such that ([F(x) — 1^(0)]^) oc jxp'*, with the Hurst exponent h > 0, was studied in 
Refs. IHinnilll] (the case h = 1/2 corresponds to the Brownian case, partly studied in Ref. [12]). 
This will not be the question of interest in the present article. Another reason for the divergence 
J dx (y{x)V{Q)) = oo is for a potential with short range correlations and large local fluctuations. 
This is the case on which we will focus here. 

For simplicity we consider a random potential V{x) uncorrelated at different positions (van- 
ishing correlation length). A model that realizes these conditions is the following random poten- 
tial 

V{x) = ^^Vn5{x - Xn) , (5) 

n 

where the weights Vn are chosen to be independent and identical random variables distributed 
according to a distribution with power law tail 
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Pi{v) OC — — for u — > ±00 (6) 



with fi > 0. Here w is a scale for the weights to make the argument of the tail dimensionless. 
The positions of impurities x„ are also chosen to be independent random variables uniformly 
distributed with a finite density p. When /i ^ 2 the second moment diverges (v^) = oo, as well 
as the correlation function of the potential since {V{x)V{x')) = p{v'^) S{x — x'). 

Some exact results have been obtained for a tight binding Hamiltonian with random on-site 
potential (Anderson model) distributed according to a Cauchy law [TSKHKlS] (§10.3 of Ref. [6] 
or Ref. [7]). Note also that fluctuations of the envelope of the wave function and conductance 
statistics for this discrete model were studied for power law disorder in the recent works [TGIITT] . 
The Anderson model (AM) can be mapped [7j onto the problem we are interested in here for 
fixed impurity positions and with /i = 1. Potentials for fixed and random impurity positions 
share some features, however, when impurities of random weights form a lattice, x„ = n/p, the 
trace of the lattice remains for arbitrary large energies (band edges remain at /c„ = mrp). This 
makes the definition of a high energy regime less convenient. 

Ricatti variable.— We follow the ideas introduced in Ref. [18] in order to study the spectrum, 
and apply them to the localization problem. Let us introduce the Ricatti variable z{x) = ^^^-j ■ 
From the Schrodinger equation we see that it obeys a Langevin equation z' = —E — + V{x) 
for initial condition z{{}) = oo. The distribution T(z; x) of the Ricatti variable obeys the integro- 
differential equation 



^T{z- x) = ^ [{E + z^)T{z; x)] + p j dvp{v) [T{z -v'x)- T{z; x)] . (7) 



The first term in the right hand side is the drift term coming from the force —{E -\- z'^) and 
the second a jump term originating from the random potential ([5]). For x sufficiently large, 
the distribution reaches a limiting distribution T(z) for a steady current [6l[7]. Current of the 



the normalized wave functions ^p{x) of the stationary Schrodinger equation on a finite interval [0, L] satisfy- 
ing boundary conditions v?(0) — ip{L) = 0, what can only be achieved for a discrete set of energies (Sturm- 
Liouville problem). From this scheme we expect that the normalized wave functions present the structure 
~ sin(A;j: + So) 6"'^"^"'^^'°'=. Note however that this simple picture neglects the important fact that, in 
the exponential, ^{x) has large absolute fluctuations despite it presents negligible fluctuations relatively to its 
average when J dx (V^(a;)V^(0)) < oo. These fluctuations play a very important role since they induce large 
fluctuations of the normalization of the wave function (see §13.3 of Ref. [6]) [8]. 
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Ricatti variable through M gives the number of zeros of the wave function per unit length. This 
is also the integrated density of states (IDoS) per unit length N(E), therefore 



N{E) = {E + z^)T{z) - p / dvp{v) / dz' T{z') 



(8) 



Imposing normalization of the solution of this integral equation gives the IDoS. Knowing the 
limiting distribution T{z), the Lyapunov exponent can be obtained from [6] 7 = (z). Since 
T(z — > ±00) ~ N(E)/z'^, in order to deal with well defined integral it is understood that 
calculation of the Lyapunov exponent involves the antisymmetric part of the distribution : 7 = 
Jdzzl[T{z)-T{-z)]). 

Let us study the high energy Lyapunov exponent. For that purpose we solve the integral 
equation ([8]) by perturbation starting from the solution for V{x) = 0. In the absence of disorder 



{p{v) = 6{v)) we have To{z) 



We expand the distribution T{z) = To{z) + Ti{z) + 



in powers of the density p, as well as the IDoS. Then eq. ([8]) is solved recursively order by order. 
We easily obtain Ti{z) from which we deduce 



7(£; = A;^ ^ 00) 
This gives the general formula 



dz- 



dvp{v) 



arctan ■ 



arctan ■ 



In 
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(9) 



(10) 



where the averaging is now taken over the 6-peak weights f„. This is the first term of a "concen- 
tration expansion" that can be systematically developed [6J (eq. (|1U|) was derived in the §10.4 
of this latter reference for non random weights Vn. Additional averaging in eq. (jlOp follows from 
the property of additivity of the variable (,{x)). 



New energy dependences.— We first consider the high energy Lyapunov exponent, \/E = 
k p, w, when the weights are distributed according to ([6]). We write pi{v) = ^f{v/w) 
where /(y) is a dimensionless symmetric function such that f{y — > ±00) ~ C\y\~'^~'^. We 
divide the integral 7 ~ p dy f{y) ln[l + (^y)^] into three parts : 7 ~ p[(^)^ <iyf{y) v'^ + 



■■2k 



fC /f dy y^-^ + 2C dy y-^-^ Hf,y 



2k /w 



• For p > 2 the Lyapunov exponent is dominated by smallest y (< 1). We obtain 7 oc p{^)'^ 
that corresponds to expand the logarithm of eq. (jlOp for small v. This is the result of eq. @ : 

• For p = 2, eq. ^ cannot be used since (u^) = 00. The integral giving the Lyapunov exponent 
is dominated by the intermediate scale 1 ^ y ^ k/w : 



For < /i < 2, the Lyapunov exponent is dominated by largest y (> k/w). We obtain 



.2^ fwy 



(11) 



(12) 



The numerical dimensionless prefactors depend on the precise form of the distribution and not 
only on its tail. 

The fluctuations of the random weights can be further increased by considering a distribution 
with tail : 



P2{v) OC 



1 

In^+'^l^l 



for V ±00 



(13) 
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for /i > 0. When /i > 1 we find that the Lyapunov exponent decays logarithmically with energy 



7(A;2) . (14) 

The case < ;U ^ 1 is discussed in the next section. 

Super localization.— On the other hand, for distribution (|13p with ^ 1, not only the second 
moment diverges (w^) = oo, but the expression ([TO]) shows that the Lyapunov exponent diverges 
as well : 7 = 00. This indicates that the logarithm of the envelope of the wave function, (,{x), 
presents different scaling properties with x. In order to analyze this, we remark that the variable 
^(x) is constant between two impurities and makes a jump = S,{x^) — ^(x^) ~ ln|f,i| 
across the impurity n (see below, the section on numerics). Therefore .^(x) behaves as the 
sum of N px independent variables, each distributed according to a power law distribution 
p(A^) oc A^""'^"'^. Using well known results (recalled in appendix A) we obtain 

C(x) ~ {px)^/^' for < /i < 1 (15) 

~ (px) In(px) for /X = 1 . (16) 

The envelope of the wave function presents a decay faster than a simple exponential. This 
phenomenon is called superlocalization and has been recently studied for a discrete model in 
Ref. |l9pl . Characterization of the localization properties cannot be limited to the typical be- 
haviours (|15|16p since the variable £,{x) presents large fluctuations. Its distribution is character- 
ized by a power law tail 

V{^;x)^l/e'-^ (17) 

with the same exponent as the one involved in the distribution of the weights. 

Conductance. - We give another interpretation of the previous result in terms of the conductance 
of a finite disordered interval of length L. The dimensionless conductance is equal to the trans- 
mission probability (Landauer formula), and presents the same exponential decay as the square 
of the wave function modulus. Therefore we can write g ~ e~^^^^^ (a more precise definition of 
the reflection coefficient within the phase formalism can be found in Ref. [5]). 

Let us first recall some well known results valid for a potential with finite local fluctuations. 
At high energy, when dH) holds, ^{x) behaves like a Brownian motion with driflU [5] : ^(x) 
7 x+^W{x), where W{x) is a Wiener procesfl It follows that the distribution of the logarithm 
of the conductance is Gaussian n(ln g) ~ ^q^^l ^^P ~s^i^^9~'^'^^^)'^ [^^■ The typical value of 
the conductance is gtyp ~ e"^"^^ (while (ln5()typ ~ — 27L) however fluctuations of the logarithm 
are associated to a much larger scale gguct ~ e~'^^^. 

Distribution of the conductance in the Anderson model with power law disorder has been 
studied in Ref. [17] where some power law distribution of the conductance was obtained for g ^ 0. 

In the superlocalization regime, the behaviour (jlSp is associated with a decay of the conduc- 
tance g ~ e~^^^'^ . The distribution (I17p can be related to the conductance distribution : 

n(lnc/) ~ ,^ ^, ^. (18) 

for < /i ^ 1. 



^ Note also that such superlocahzation ^{x) ^ x^^^^'^ occurs for self affine random potentials characterized by 
long-range correlations ([V{x) - V(Q)f) oc |a;p'' with h>Q [IT] . 

^The fact that drift and the Gaussian fluctuations involve the same parameter is refered to as '^single parameter 
scaling" [20] (see also [16]) ; it holds only at high energy since it relies on the decoupling between a fast variable 
(the phase 0(x) introduced above) and the slow variable £,{x). 
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a normalized free Brownian motion such that {W{x)) = and {W{x)W{x')) — min(a::, 
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Localization for supersymmetric Hamiltonian. We consider another class of random 
Hamiltonians with the so-called supersymmetric structure 

d^ 

H = --^ + (Pixf + ^'{x) with(P{x) = ^r]nS{x-Xn) , (19) 

n 

where rjn are dimensionless uncorrelated weights, each distributed according to a distribution 
p{r]). This Hamiltonian is interesting since it presents rather different spectral and localization 
properties (in particular it leads to a delocalization transition as E ^ 0). It is related to several 
other problems as well. For example it is the square of a Dirac Hamiltonian with a random 
mass (j){x), introduced in various contexts of condensed matter physics ; the problem can also 
be related to classical diffusion in a random force field (see Refs. [HI [22l [23] for a review). We 
can follow the same strategy : the Ricatti variable z = ^ — cp obeys the Langevin type equation 
z' = —E — — 2z(/)(x) with multiplicative noise. Limiting distribution of the Ricatti variable 
for a steady current —N(E) obeys the integral equation 

dr]p{r]) J dz'T{z') . (20) 

The Lyapunov exponent is now given by 7 = {z) + {(p) . Following the same perturbative approach 
as before, we obtain 

7(£' ^ 00) ~ p (Incoshr?)^ . (21) 

For the supersymmetric Hamiltonian, when the average exists, the Lyapunov exponent reaches 
a finite value at high energy (in contrast with the decrease of the Lyapunov exponent for the 
Schrodinger Hamiltonian) : for |r/„| <C 1 it takes the form ^{E 00) ~ i J dx {(j){x)(j){0)) . 

Let us consider a power law distribution of weights p{r]) oc l/|r/|^"*''^. We see from (12 ip 
that we have 7 = 00 in this case for // ^ 1. The reason is similar to the one discussed in 
the previous paragraph for the Schrodinger equation with weights distributed according to p3|) . 
Here the variable .^(x) jumps by A,^„ ~ \r]n\ across the impurity (see below). Therefore, for 
the supersymmetric case, a power law distribution of the weights leads to the superlocalization, 
eqs. (I15I16I17P . 

Numerical calculations.— We can easily study the evolution of the phase and envelope vari- 
ables (j2|3p numerically. We denote by 9^ = 0{xn) and = ^i^n) the value of the phase 
and the envelope just before and right after the n-th (5-peak. Between two impurities we have 
0~^i — On = kin and £,^+1 ~ = 0- The length in = Xn+i — Xn denotes the distance between 
consecutive impurities. It is distributed according to a Poisson law p{i) = pe~P^. The evolution 
of the random variables across an impurity depends on the form of the random potential. We 
introduce the notation = ~ ?n = ^n+i ~ ■ 

• For the Schrodinger Hamiltonian with potential ([5]) phase evolution is given by cotg 9n — 
cotg6';;; = ^ and evolution of the envelope by = In = i ln[l + ^ sin26'~ + ^ sin^ 0-]. 

• For the supersymmetric Hamiltonian (fT9|l we have tan 6^ = e^''" tan 6" and = | l^i = 
i In [e^'?" sin2 - + e- ^'(n cos2 ] . 

IDoS is given by N{E) = limi_>oo and the Lyapunov exponent by "y{E) = lim^^oo 
To be precise we consider a specific distribution with power law tail : 





— 1 

w 1 




TTW (1 + 


1 ^ 1 
1 w 1 





This choice has the advantage that it is very easy to simulate since the cumulative distribution 
is straightforwardly obtained. Using eq. (llOp . we get the high energy Lyapunov exponent : 
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7 — sin(^^) ^(^)^ for < < 2 and 7 ~ f ) for /i = 2 (both expressions are valid for 
k ^ p, w). The case n = 1 corresponds to a Cauchy law. The right hand side of eq. (jlOp can 
be computed easily in this case and we obtain the expression 7 ~ pln{l + ^) valid in a broader 
range of energy (for k ^ p but w arbitrary) ; at high energy we recover the known energy 
dependence 7 ~ p^ (it can be recovered from discrete models |14t [7]). These expressions are 
compared to the numerical results on figured] and work perfectly well. 
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Figure 1: For the Schrddinger Hamiltonian : Lyapunov exponent as a function of the energy for 
weights Vn distributed according to 12^) . Numerical results (continuous lines) are compared with 
high energy expressions (dashed lines) derived in the text (no fit) for p = 0.34, 1 and 2. Other 
parameters are w = 1, p = 0.01 and number of impurities N = 10^. Inset : E^{E) is plotted in 
semilog scale for p = 2 in order to check its logarithmic behaviour (dashed line corresponds to 
eq. m)- 

Next we analyze the superlocalization regime : we consider the supersymmetric Hamiltonian 
for weights r]n distributed according to a law similar to (p2]) for < ^ < 1. The distribution 
■p(^; x) is plotted for different values of x on figure [2j In the inset the axes are rescaled in order 
to check that, according to (llSp . the distribution has the form 

where w^C) is a dimensionless function. After rescaling we see that the four curves corresponding 
to different values of x perfectly collapse onto each other, apart for small deviations correspond- 
ing to the smallest values of ^ and x. Finally we check that the tail of the distribution is indeed 
a power law, eq. (fTTj) : in the inset of figure [2] rescaled distributions are plotted on a log-log 
scale with w{(^) oc (^~^~'^. The agreement seems excellent. 

Conclusion.— We have analyzed the high energy localization length for random potentials with 
short range correlations and large local fluctuations such that the well-known result leading 
to £ioc oc E, is not valid. We have studied localization for potentials made of superposition 
of (5-peaks. Performing a concentration expansion, we have obtained two general high energy 

2 

formulae for the Lyapunov exponent : 7(-E) — p(ln[l-|- [6j for the Schrodinger Hamiltonian 
and ^{E) ~ p(lncosh7/)^ for the supersymmetric Hamiltonian. These formulae have been used 
to analyze the case of potential with large local fluctuations. 

For the Schrodinger case, we have shown the relation between the distribution of the weights 
of the (5-peaks and the energy decay of the Lyapunov exponent (inverse localization length £ioc)- 
Sufficiently large fluctuations of the weights, such that (f^) = 00, lead to a stronger localization 
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Figure 2: Superlocalization for the supersymmetric Hamiltonian. Distribution of the variable 
^(x) for different values of x = L/4, L/2, 3L/4 and L. Parameters are : k = 10, p = 1, A = 0.1 
(typical scale for weights rjn) and /i = 0.5. Number of impurities is N = 10^ and L = 1000. In 



the inset, straight line corresponds to eq. (17) 



effect characterized by an increase of i\oc with energy slower than hnear. These results are 
summarized in tabled) 



Potential distribution Localization length 



{v'i) 


< oo 






Aoc « E [Q\ 


p{v 


±Q0) 


oc 1/1^1'^+^ 


fi = 2 
< /!< 2 


iioc (xE/lnE 
iioc oc E^^/^ 


p{v 


±oo) 


oc ^In-^-'^ 1^1 

\v\ \ w I 


/i > 1 


Aoc oc In'^-^^; 






superlocalization 


r ^ = 1 

\ < /i < 1 


^(x) ~ xlnx 



Table 1: Energy dependence of the localization length for the Schrddinger Hamiltonian with 
random potential V{x) = X^„'Un(5(x — x„) for different broad distributions of the weights Vn- 

The understanding of fluctuations of the variable .^(x) (i.e. of the localization length) plays 
a major role to analyze universal statistical properties of Wigner time delay [8]. It would 
be an interesting issue to study how the statistics of Wigner time delay are affected by the 
unconventional localization properties analyzed here. 

Acknowlegments.— We thank Alain Comtet and Satya Majumdar for interesting remarks and 
Jean-Marc Luck for bringing to our attention Refs. |11[ I19j . 

Appendix A.— We recall well known results on the distribution of the sum of independent 
and identically distributed (i.i.d) random variables. Let us consider N i.i.d. positive variables 
y„ and their sum = Yln=i V^- If {Un) < oo the statistical properties of Yat are given by 
central limit theorem for N ^ oo : Gaussian distribution centered on {Yjy) = N{y) of variance 
(Y^)^ = (Y^) - (yjv)2 = N{y^)c. If the distribution of the y„'s presents a power law tail 
p{y) oc l/y^+^ with < ;U ^ 2 such that (y^) = oo, the situation is different : 
• For < ;U < 1 all moments of yn diverge. Let us consider the characteristic function g{p) = 
(^Q-pyy \yg write g{p) = I — dy (1 — e~^^)p(y) from which we see that g{p — > 0) ~ 1 — Cp^ 
where C is some constant related to the prefactor of the power law tail olp{y). The characteristic 
function for Y/v reads Gn{p ^ 0) ~ e~^^^'^ . This shows that the related distribution P/v(y) 
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presents a similar power law tail and involves the typical scale Yiv ~ N^^'^- 

• For = 1. A similar analysis gives g{p ^ 0) ~ 1 — Cpln 1/p and therefore ~ A^ln A^. 

• For 1 < < 2 the first moment is finite (Itv) = N {y) however fluctuations are larger than in 
the normal case {Y^)c ~ iV^/^. 

• For fi = 2 fluctuations are {Y^)c ^ NlnN. 

• For n> 2, central limit theorem applies. 
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